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Q\ . Medium modification of jet fragmentation functions and parton energy loss in cold and 
^ ! hot matter are reviewed. The predicted nuclear modification of the jet fragmentation 
>- \ function agrees well with the recent HERMES data with a resultant energy loss dE/dx m 
0.5 GeV/fm. From the the recent PHENIX data of high ttq spectra in central Au + Au 
collisions at -/s = 130 GeV, one extracts an energy loss for a 10 GeV parton that is 
equivalent to dE/dx = 7.3 GeV/fm in a static medium with the same gluon density as 
! in the initial stage of the collision at r = 0.2 fm/c. Constraints on jet quenching by the 
central rapidity density of charged hadrons is also discussed. 
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1. INTRODUCTION 



a: 

It is with a great respect that I present this talk on this special occasion of celebrating 
Helmut Satz's scientific career and his unwavering role in high-energy heavy-ion physics. 
I first met Helmut when I was still a graduate student at my first Quark Matter meeting 
in 1988. His name was already synonymous to quark- gluon plasma and J/^f suppression 
for me at that time. Because of our common interests in hard processes in heavy-ion 
collisions, we started working together in 1994 to coordinate a working group called Hard 
Probes Collaboration. Since then this working group has become a forum where we meet 
regularly together with other colleagues to discuss problems related to hard processes 
in high-energy heavy-ion collisions. Helmut has been an inspiration for all of us in the 
working group and in the community of high-energy nuclear physics. 

The subject of jet quenching and parton energy loss is one aspect of hard processes 
in heavy-ion collisions. However, with the advent of the RHIC facility and experiments, 
jets of high-energy partons will become an important and useful tool to the study of the 
properties of dense matter formed in high-energy heavy-ion collisions. Because large pt 
partons are produced very early in heavy-ion collisions and their production rates can 
be calibrated in pp and pA collisions at the same energy, they are ideal probes of the 
dense matter that is formed in the same reaction. What probes the dense medium is the 
scattering induced energy loss suffered by an energetic parton as it propagates through the 
matter. The parton energy loss is directly related to the parton density of the medium. 

Theoretical studies of the parton energy loss in hot medium date back to the first 
attempt by Bjorken [|TJ to calculate elastic energy loss of a parton via elastic scattering in 
the hot medium. A simple estimate can be given by the thermal averaged energy transfer 
v e \ ~ <?i/2w of the jet parton to a thermal parton with energy lu, q± being the transverse 
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momentum transfer of the elastic scattering. The resultant elastic energy loss 

dx 2 \2^ ) K ' 

is sensitive to the temperature of the thermal medium but is in general small compared 
to radiative energy loss. Here, \x is the Debye screening mass and C<i is the Casimir of 
the propagating parton in its fundamental presentation. The elastic energy loss can also 
be calculated within finite temperature QCD || with a similar result, but with a more 
careful and consistent treatment of screening effect. 

Though there had been estimates of the radiative parton energy loss using the uncer- 
tainty principle a first theoretical study of QCD radiative parton energy loss incor- 
porating Landau-Pomeranchuk-Migdal interference effect || is by Gyulassy and myself 
[[J where multiple parton scattering is modeled by a screened Coulomb potential model. 
Baier et al. (BDMPS) later considered the effect of gluon rescattering which turned 
out to be very important for gluon radiation induced by multiple scattering in a dense 
medium. These two studies have ushered in many recent works on the subject, including 
a path integral approach to the problem || and opacity expansion framework PJITfl which 
is more suitable for multiple parton scattering in a thin plasma. The radiative parton 
energy loss to the leading order of the opacity n = L/X in the thin plasma of size L is 
estimated as |9|JTT[] 



dE iad a s /i 2 L ( 2E\ 

— ^—-In^—j, (2) 

where A is the gluon's mean-free-path in the medium. The unique L dependence of the 
parton energy loss is a consequence of the non-Abelian LMP interference effect in a QCD 
medium. It is also shown in a recent study |TT| that thermal absorption and stimulated 



emission in a thermal environment can be neglected for high energy partons (E ^> fx) 
while they are important for intermediate energy partons. 

Using this latest result one can estimate the total energy loss for a parton with initial 
energy E = 40 GeV to be about AE 10 GeV after it propagates a distance of L = 6 fm 
in a medium with /i = 0.5 GeV and A = 1 fm. For an expanding system, the total energy 
loss is reduced by a factor of 2t /L from the static value |T2"| , |2"T| 1 . Assuming that most of 



this energy loss is carried by gluons outside the jet cone |T3, measuring the energy loss 



would require the experimental resolution SE to be much smaller than the total energy 



loss AE. With the measured total multiplicity density dN/drj « 900 |H| and energy 



density dEr/dr/ m 500 GeV |16 in central Au + Au collisions at y/s = 130 GeV, one 
can estimate that the average total background energy within the jet cone (5rj = 1 and 
5(f) = 1) is about T,E T m 80 GeV with a fluctuation of 5E T « 10 GeV. It is therefore 
very difficult, if not impossible, to determine the energy of a jet on a event-by-event base 



17j • Since high px hadrons in hadron and nuclear collisions come from fragmentation of 
high Et jets, energy loss naturally leads to suppression of high px hadron spectra. Miklos 
Gyulassy and I then proposed [IIJ that one has to reply on measuring the suppression of 
high pt hadrons to study parton energy loss in heavy-ion collisions. Since inclusive hadron 
spectra is a convolution of jet production cross section and the jet fragmentation function 
in pQCD, the suppression of inclusive high p T hadron spectra is a direct consequence of 
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the medium modification of the jet fragmentation function induced by parton energy loss. 
Assuming that jet fragmentation function is the same for the final leading parton with a 



reduced energy, the modified fragmentation function can be assumed as ]T9] 



Therefore, in this effective model the measured modification of fragmentation function 
can be directly related to the parton energy loss. 

2. MODIFIED FRAGMENTATION FUNCTIONS 

Since a jet parton is always produced via a hard process involving a large momentum 
scale, it should also have final state radiation with and without rescattering leading to 
the DGLAP evolution equation of fragmentation functions. Such final state radiation 
effectively acts as a self-quenching mechanism softening the leading parton momentum 
distribution. This process is quite similar to the induced gluon radiation and the two 
should have strong interference effect @,|nj. It is therefore natural to study jet quenching 
and modified fragmentation function in the framework of modified DGLAP evolution 
equations in a medium [ 2(| . 



The simplest case for jet quenching is deeply inelastic scattering of an electron off a nu- 
cleus target where the virtual photon knocks one quark out of a nucleon inside a nucleus. 
The quark then will have to propagate through the rest of the nucleus and possibly scat- 
ter again with other nucleons with induced gluon radiation. The induced gluon radiation 
reduces the quark's energy before it fragments into hadrons with a modified fragmenta- 
tion function. One can study the nuclear modification of the fragmentation function by 
comparing it with the same measurement in DIS with a nucleon target. 

We work in an infinite momentum frame, where the photon carries momentum q = 
[— XbP + , q~ , 0_i_] and the momentum of the target per nucleon is p = [p + ,0,0j_] with the 
Bjorken variable defined as xb = Q 2 /2q~p + . The differential semi-inclusive hadronic 
tensor in a collinear factorization approximation to the leading twist can be written as 

dW s 

dz h 

where f^(x,fij) is the quark distribution of the nucleus, H^ v {x,p,q) is the hard part 
of 7* + q scattering and D q _>h(zh, fi 2 ) is the quark fragmentation function in vacuum as 
measured in the e + e~ annihilation process. The scale ft 2 dependence of the fragmentation 
function comes from the radiative correction as shown in Fig. |l| Taking into account of 
both the radiative (central cut diagram) and the virtual correction (virtual cut diagram), 
the renormalized fragmentation function is defined as, 

, o \ j-. / % dtt a s f 1 dz 

D q ^ h {Z hl fl ) EE D q ^ h (z h ) + / ^2-—/ — 

X bq~,qg{z)D q ^ h (z h /z) +lq^gq(z)D g ^ h (z h /z)} , (5) 

which satisfies the DGLAP evolution equation, 
dD q _+ h (z h , fi 2 ) _ a s f 1 dz 



J2 J dxf^(x, fi 2 )H^(x,p,q)D q ^ h (z h , fi 2 ) , (4) 



d In /i 2 27r 



l q -> qg (z)D q ^ h (z h /z } fi ) +7 q ^ gq (z)Dg^ h (z h /z,n ) .(6) 
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Here ^ qg (z) = 7™(1 - z) = C F [{1 + z 2 )/(l - z) + + (3/2)5(1 - z)\ is the splitting 
function. 

In a nucleus target, the outgoing quark can scatter again with another parton from the 
nucleus. The additional scattering may induce additional gluon radiation and cause the 
leading quark to lose energy. Such induced gluon radiation will effectively give rise to 
additional terms in the evolution equation leading to modification of the fragmentation 
functions in a medium. Contributions from multiple parton scattering are always non- 
leading twist. However we will consider only those that are enhanced by the nuclear size 




Figure 1. Diagram of final state radiation 
correction to the jet fragmentation function 
with central and virtual cuts. 



Figure 2. One example diagram of gluon 
radiation associated with double scatter- 
ing with different cuts (central, left, right 
and virtual). 



We work in the LQS framework jn|] in which the twist-four contributions can be ex- 
pressed as the convolution of partonic hard parts and four-parton matrix elements. At 
the lowest order (processes without gluon radiation) in this framework, rescattering with 
collinear gluons gives the eikonal contribution to the gauge-invariant leading-twist and 
lowest-order result, assuming collinear factorization of the quark fragmentation function. 
The radiative correction to double scattering processes with another gluon from the nu- 
cleus, shown in Fig. [2] as an example, generally involves matrix elements of four parton 
fields and the cross section of four parton scattering. In the LQS framework, one makes 
collinear expansion of the four-parton scattering cross section in terms of the transverse 
momentum of the scattering gluons. The first term in the collinear expansion corresponds 
to rescattering with collinear gluons. Similar to the leading order case, it contributes only 
to the eikonal correction of the radiative processes at the leading twist as a consequence of 
the cancellation between different cut diagrams (central, right and left cut). The leading 
twist-four contribution then comes from the second derivative of the four-parton scat- 
tering cross section in the central-cut diagrams. With the quadratic term k^_k±, this 
contribution depends on the gluon field strength, because 
partial integrations. 

Considering all possible diagrams in addition to the example shown in Fig. ||, including 
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the virtual corrections (virtual cut), one get the leading twist-four contribution to the 
semi-inclusive cross section from double scattering, 



W D 

dzh 



d^_ 



1 dza s 
Zh z 2n 



x 



D n 



+ D Q 



where 



*h{Zh/z) 

,{zh/z) 
dy 



1+* r-nA 



1 - z. , 
l + (l-z) 2 

z+ 



T*{x,x L ) + 5(z-l)M±{x,%) 



T A (x J x L ) + 6(z)AT A (xJ 2 T 



f JL- dy^dy^e i{x+XL)p+y +iXTP fo i } (1 - e ~ iXLP+ ^ )(1 - e 
J 2ir 



(7) 



-ix L p+(y )\ 



2tt 

(A|^(0) 7 



^ + (y 2 *)^(yr)V' < ,(y-)l^>e(-yj)e(y--y 1 - 



is the quark-gluon correlation function which essentially contains four independent twist- 
four parton matrix elements in a nucleus. Here xl = £ 2 ± /2p + q~ z(l — z) and xt = 
(kj) /2p + q~ z(l — z). In the central-cut diagrams (Fig. |2] for example) where the lead- 
ing higher-twist result comes from, there are typically four contributions from each cut 
diagram. The radiated gluon can either come as the final state radiation of the 7-quark 
scattering or as the initial state radiation of the secondary quark-gluon scattering. The 
amplitudes of these initial and final state radiation have opposite signs with also a phase 
difference xip + y~ ■ The sum of these two radiation processes and their interferences gives 



rise to the dipole-like factor (1 — e" 



ix L p+y 2 



)(1 



-ix L p+(y -y x )\ 



in the effective two-parton 



correlation function that enters into the double-scattering cross section with induced gluon 
radiation. This is exactly the LPM effect of bremsstrahlung in medium. 

Summing up the single and double scattering contributions to the semi-inclusive pro- 
cess, one can define the effective modified fragmentation function as 



dW, 



dz h 



dW s 

""'IMS 

dzh. 



+ 



dW, 



D 
I"' 



dzh 



E / dx f^( x ^ 2 i) H ^( x ^P^) D g-^h(z h ,fi 2 ), 



(9) 



(8) 



D q ^ h (z h ,n ) 
AD q ^ h (z h , fx 2 ) 



D q ^ h (z h , fi 2 ) + AD q ^ h (z h , /i 2 ) 



f di 2 T a s r 1 dz 



l$r 2tt Jz h z 



A^ qg (z)D q ^ h A + A^ gq (z)D g ^ h {-^' 



where D a ^ h (zh, H 2 ) is the leading-twist parton fragmentation function and 



A^ qg (z) 
A lq ^ gq (z) 



1 + z 2 



L-Z) + 

A 7 g_ 99 (l - Z). 



T A (x,x L ) + 6(l-z)AT q A (x,e 2 T 



C a 2-kol, 



{e T +{k 2 T ))N c fMx^ 



:i2) 



are the modified splitting functions for the induced gluon radiation. Similar as the vacuum 
case, the ^-function part is from the virtual correction contribution with AT A (x, i\) 
defined as 

,.1 1 

(13) 



AT A (x,£ 2 T ) = J dz—[2^(x,x L )\ g=1 -(l + g ( l )T£ l (x,x L ) 
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Such virtual corrections are important to ensure the infrared safety of the modified frag- 
mentation function and the unitarity of the gluon radiation processes. This virtual cor- 



rection is equivalent in nature to the absorption processes in some effective models [22 



To evaluate the modified fragmentation, one needs to know the two-parton correlation 
function T a (x,xl) which consists of both diagonal and off-diagonal twist-four parton 
matrices. We generalize the factorization assumption of LQS |^TJ to both types of matri- 
ces. Assuming a Gaussian nuclear distribution in the rest frame, p(r) ~ exp(— r 2 /2R A ), 
R A = 1.12 A 1 / 3 fm, we express T A in terms of single parton distributions, 

T q A g (x,x L ) = Cm N R A f q A (x)(l - e-*M), (14) 

where xa = 1/ttinRa, and is the nucleon's mass. The off-diagonal terms involve 
transferring momentum xl between different nucleons inside a nucleus and thus should be 
suppressed for large nuclear size or large momentum fraction x^. Notice that Tf = l/xip + 
is the gluon's formation time. Thus, xl/xa = La/tj with La = Ra^nIv^ being the 
nuclear size in the infinite momentum frame. 

Because of the LPM interference effect, the above effective parton correlation and the 
induced gluon emission vanishes when xl/xa <C 1. Therefore, the formation time of 
the gluon radiation due to the LPM interference requires the radiated gluon to have a 
minimum transverse momentum ij, ~ Q 2 /MRa ~ Q 2 /A 1 / 3 . The nuclear corrections 
to the fragmentation function due to double parton scattering will then be in the order 
of asA 1 / 3 /£j, ~ a. s A 2 l 3 /Q 2 , which depends quadratically on the nuclear size. For large 
values of A and Q 2 , these corrections are leading and yet the requirement i 2 ^ <C Q 2 for the 
logarithmic approximation in deriving the modified fragmentation function is still valid. 




Figure 3. The predicted nuclear modification 
of jet fragmentation function is compared to 



the HERMES data 23 



Figure 4. Comparison of the calculated 
nuclear modification with the effective 
model in Eq. @ with Az = 0.6(z g ). 



Shown in Fig. 3 are the calculated nuclear modification factor of the fragmentation 
function for 14 N and SA Kr targets as compared to the recent HERMES data [ 23| . The 
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predicted shape of the z dependence and the quadratic nuclear size dependence agrees 
well with the experimental data. The energy dependence of the suppression also has 
excellent agreement with our prediction [p4 |. What is amazing is the clear quadratic A 2 ^ 3 



nuclear size dependence of the suppression which is a true QCD non-Abelian effect. In 
fitting the data of the overall suppression for 14 N target we obtain the only parameter 
in our calculation, Ca 2 = 0.00065 GeV 2 . This parameter is shown to be related to 



the transverse momentum nuclear broadening of Drell-Yan dilepton in pA collisions p5 
(Aq 2 ± ) = Cna s /N c x A . With an experimental value f25[ of (Aqj) = 0.016A 1 / 3 GeV 2 and 



with a s = 0.21 (at M 2 = 40 GeV 2 ), one finds Ca 2 = 0.00057 GeV 2 . This indicates that 
C might have some Q 2 dependence. With this value of C we can also predict the nuclear 
transverse momentum broadening of single jets in DIS p3. 

With the modified fragmentation function in Eq. (|10D , one can calculate theoretically 
the average energy loss by the quark, which is the energy carried away by the radiated 
gluons, 

AE = v(Az g } » Ca 2 m N R 2 A (C A /N c )3 ln(l/2s B ). (15) 



With the value of a 2 C, and L A = R A V2n one gets the quark energy loss dE/dx ~ 0.5 
GeV/fm for a Au nucleartarget. 

3. JET QUENCHING IN HEAVY-ION COLLISIONS 

In high-energy heavy-ion collisions, the jet production rate is not affected by the for- 
mation of dense matter and the final state multiple scattering. One can assume that the 
high px hadron spectra can then be given by the convolution of the jet production cross 
section and the medium modified jet fragmentation function D h / C (z c , Q 2 ), 



do h AB 
dyd 2 pT 



K J d 2 bd 2 rt A (r)t B (\h - r\)Y^ J dx a dx b d 2 k aT d 2 k bT g B (k bT , Q 2 ,\h 



abed 

x gA(k aT , Q 2 , r)f a/A (x a , Q 2 , r)f b/B (x b , Q 2 , |b - Y \) Dh / c<yZc ' ® ) ^1^1 ,6) 

TIZ C dt 



where t A (r) is the thickness function of the nucleus A, f a /A(x a ,Q 2 ,r) is the parton dis- 
tribution in a nucleus, g A (k a T,Q 2 ,r) is the distribution of parton intrinsic transverse 
momentum with nuclear broadening ||26|| . 

In principle, one should use the modified fragmentation function evaluated according 
to the pQCQ calculation for a dense medium. However, before that can be done in a 
practical manner, we have used the effective approach in Eq. (|3D by rescaling the fractional 
momentum by 1—Az to take into account of the parton energy loss. To verify whether such 
an effective approach is adequate, we compare the two modified fragmentation functions 
in Fig. |1 We found that the effective model (dashed lines) can reproduce the pQCD 
result (solid lines) of Eq. (|I0| ) very well, but only when Az is set to be Az ps 0.6(z g ). 
Therefore the actual averaged parton energy loss should be about 1.6 times of that used in 
the effective modified fragmentation function. This difference is caused by the absorptive 
processes or unitarity correction effect in the full pQCD calculation. 

Unlike in DIS nuclear scattering, the dense medium in high-energy heavy-ion collisions 
is not static. It has to go through rapid expansion which should also affect the effective 
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total parton energy loss. The total energy loss extracted from experiments should be a 
quantity that is averaged over the whole evolution history of the expanding system. It is 
therefore useful to convert the averaged quantity to an energy loss in a static system that 
has the same parton density as the expanding system at its initial stage. If the averaged 
total parton energy loss in a longitudinally expanding system with a transverse size R is 
AEia, one finds |Z7[] that the corresponding parton energy loss in a static system with the 
same initial parton density would be AE = AEi^(R/2tq). Here r is the initial formation 
time of the dense medium. 



la □ UA1 Vs=200 GeV 



p + pbar -> (h + + h")/2 : 




Pt (GeV/c) 




Vs (GeV) 



Figure 5. pQCD parton model calculation 
of the charged hadron and pion spectra in pp 
and central Au+Au collisions compared with 
the experimental data [2S,2?J. The effective 
modified fragmentation function is used in 
the calculation. 



Figure 6. The energy dependence of the 
charged multiplicity density per partici- 
pant nucleon pair in the two-component 
model with and without parton shadow- 
ing compared with the experimental data 

mm 



Comparing the recent PHENIX data with Eq. fll6|) as shown in Fig. || one can 
extract a value of dE/dx = 0.23 \aE/ ln5GeV/fm with a logarithmic energy dependence 
that one needs to use in the effective modified fragmentation function in fitting the data. 
Taking into account the unitarity correction effect and the expansion, this corresponds 
to an effective energy loss dE/dx = 1.6 x 0.23(-R/2to) In E 1 / In 5. in a static system with 
a density similar to the initial stage of the expanding system at To. With R ~ 6 fm and 
To ~ 1/po — 0.2 fm, this would give dE/dx ~ 7.3 GeV/fm for a 10 GeV parton, which 
is about 15 times of that in a cold nuclear matter as extracted from the DIS processes. 
Since the parton energy loss is directly proportional to gluon density, this implies that 
the gluon density in the initial stage of Au + Au collisions is about 15 times higher than 
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that inside a cold nucleus. 

4. GLOBAL CONSTRAINTS OF PARTON ENERGY LOSS 

When a fast parton experiences multiple scattering and induced gluon radiation, the 
radiated gluon would also contribute to the final parton production leading to the en- 
hancement of modified fragmentation function at small z as shown in Fig. f|. These 
partons might get thermalized in the medium, but will contribute to the final hadron pro- 
duction. In the HIJING model with the default setting, a simple dE/dx=2 GeV/fm 



is assumed for every jet with pt > 2 GeV/c 2 . Such a default setting predicted an en- 
hancement of the total hadron multiplicity in the central region |1| and a strong energy 
dependence [[H]] which is now excluded by the new PHOBOS data ||15[ of Au + Au colli- 



sions at y/s = 200 GeV. This either implies that the dE/dx = 2 GeV/fm used in HIJING 
is too large or the onset of jet quenching is set at too low px or both. According to the 



recent study of parton energy loss with detailed balance ||11|| , gluon absorption by the fast 
parton in the thermal medium reduces the parton energy loss significantly for low energy 
partons. Taking into account this strong energy dependence of dE/dx, one might expect 
that the threshold for gluon radiation could be higher than what HIJING used. In fact, 
if one increases the threshold to 3 GeV in HIJING, one can still fit both dN c ^/dr] and 
the suppressed hadron spectra at high p^. Since the production rate of mini-jets with 
Pt > 3 GeV/c 2 is very small, their contribution due to induced bremsstrahlung to the 
total charged multiplicity is negligible at the RHIC energies. By similar argument, the 
effect of parton thermalization on the total hadron multiplicity might also be small. 

Neglecting the effect of jet quenching on the total multiplicity, one can assume the 
final dNch/dr] to be proportional to the total number of mini-jets produced in addition 
to the soft particle production. Assuming that the mini-jet production is proportional 
to the number of binary collisions and the soft part is proportional to the number of 
participating nucleons, one get in this simple two-component model 

^ = l(N^ t )(n) s + (n),(iV binary )^^, (17) 

where a^(s) is the averaged inclusive jet cross section per NN in AA collisions. The 
energy dependence of the above multiplicity is shown in Fig. || The parameters (n) s = 
1.6 and (n)^ = 2.2 is fixed by the p + p(p) data. The average number of participant 
nucleons and number of binary collisions for given impact-parameters can be estimated 
using HIJING Monte Carlo simulation. If one assumes that the jet production cross 
section a^(s) is the same as in p + p collisions, the resultant energy dependence of the 
multiplicity density in central nuclear collisions is much stronger than the RHIC data. 
Therefore, one has to consider nuclear effects of jet production in heavy-ion collisions. 



Using a more recent parameterization ||33|| , we found that the RHIC data requires a 



stronger nuclear shadowing of gluon distributions than quarks. 
5. CONCLUSIONS 



We have calculated the medium modification of the jet fragmentation functions due to 
gluon radiation induced by the multiple parton scattering. The predictions of the shape, 
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energy dependence, and the quadratic nuclear size A 2 / 3 dependence of the modification 
agree well with the recent HERMES data. The resultant parton energy loss in the cold 
nuclear medium is estimated to be about 0.5 GeV/fm inside Kr nuclei. Comparing to 
the QCD result of the modification of fragmentation function, we found that the actual 
averaged energy loss is about 1.6 times that of the effective energy loss used in a earlier 
effective model for the same modification. Considering the effect of expansion, we found 
that the recent PHENIX data imply a medium induced energy loss in central Au + Au 
collisions equivalent to 7.3 GeV/fm in a static medium with the same gluon density as in 
the initial stage of the collision. This is about 15 times larger than the energy loss in a 
cold nucleus. Due to detailed balance in induced radiation and absorption, we argue that 
parton energy loss for small and medium high energy jets is very small. Therefore, the 
contribution from their induced radiation to the hadron multiplicity is negligible. 

This work is supported by the Director, Office of Energy Research, Office of High 
Energy and Nuclear Physics, Division of Nuclear Physics, and by the Office of Basic 
Energy Science, Division of Nuclear Science, of the U.S. Department of Energy under 
Contract No. DE-AC03-76SF00098 and in part by NSFC under project 19928511. 
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